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ABSTRACT 

A re-analysis of the Stokes-parameter spectra obtained of the z=0.692 21 cm absorption line toward 3C 286 
shows that our original claimed detection of Zeeman splitting by a line-of-sight magnetic field, Z?i os = 87 /iG is 
incorrect. Because of an insidious software error, what we reported as Stokes V is actually Stokes U : the revised 
Stokes V spectrum indicates a 3-cr upper limit of Z?i os < 17 /iG. The correct analysis reveals an absorption 
feature in fractional polarization that is offset in velocity from the Stokes / spectrum by -1.9 km s -1 . The 
polarization position-angle spectrum shows a dip that is also significantly offset from the Stokes / feature, but 
at a velocity that differs slightly from the absorption feature in fractional polarization. We model the absorption 
feature with 3 velocity components against the core-jet structure of 3C 286. Our \ 2 minimization fitting results 
in components with differing (1) ratios of H I column density to spin temperature, (2) velocity centroids, and 
(3) velocity dispersions. The change in polarization position angle with frequency implies incomplete coverage 
of the background jet source by the absorber. It also implies a spatial variation of the polarization position angle 
across the jet source, which is observed at frequencies higher than the 839.4 MHz absorption frequency. The 
multi-component structure of the gas is best understood in terms of components with spatial scales of ~ 100 
pc comprised of hundreds of low-temperature (T < 200 K) clouds with linear dimensions of << 100 pc. We 
conclude that previous attempts to model the foreground gas with a single uniform cloud are incorrect. 
Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

In order to detect magnetic fields in galaxies with signifi- 
cant redshifts, we began a search for Zeeman splitting of 21 
cm absorption lines arising in damped lyman alpha systems 
(DLAs) toward background quasars (Wolfe et al. 2005) se- 
lected to be radio bright. The detection of Zeeman splitting 
yields the strength, direction, and redshift of the magnetic 
field, which is an advantage over measurements of Faraday ro- 
tation, which yield rotation measures accompanied with huge 
uncertainties in the strength and redshift of the B fields along 
the sightlines to background quasars (Kronberg et al. 2008). 
In Wolfe et al. (2008) we reported the detection of Zeeman 
splitting in the z=0.692 DLA toward 3C 286. The purpose 
of this paper is to show that the reported detection was erro- 
neous. In $2] we describe the errors that led to the claimed 
detection and describe our new data reduction procedures. In 
<2] we correct the previous errors and give a correct descrip- 
tion of the polarization spectra. The new spectra are shown in 

m 

In the absence of detectable Zeeman splitting we can only 
set an upper limit to the magnetic field in the z=0.692 DLA 
toward 3C 286. However, in this paper we discuss how the 
spatial distribution of the strong linear polarization of 3C 286 
itself can be used in combination with observed spectral vari- 
ation of linear polarization across the absorption profile of the 
redshifted 21 cm line to infer physical properties of the ab- 
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sorbing H I gas. We demonstrate how comparison between 
properties of the linearly polarized and unpolarized absorp- 
tion spectra provides new information about the spatial struc- 
ture of this gas on scales of ~ 100 pc. We use the new data 
to determine both the kinetic temperature and hyperfine spin 
temperature of the gas, which in turn tells us about its thermal 
state. We use these results to evaluate the standard assump- 
tion of a large-scale, uniform cloud. We start by constructing 
models to describe the re-analyzed spectra in terms of multi- 
ple clouds distributed across the spatially extended structure 
of 3C 286. In $5] 1 we use VLBI maps to exhibit the core-jet 
structure of 3C 286 near the 21 cm absorption frequency. In 
$5]2 we describe a two-cloud configuration to model the ab- 
sorption feature. We use a chi square minimization technique 
to fit the model to the Stokes / spectrum and the fractional po- 
larization spectrum and show that the model does not work. 
In $5]3 we show how an adjustment in covering factors and 
the presence of velocity gradients in the cloud toward the po- 
larized jet source results in a three-component model that pro- 
vides an adequate fit to the Stokes /, fractional polarization, 
and polarization position-angle spectra. The implications of 
the results are discussed in §|6] All the results are then sum- 
marized in §0 

Throughout this paper we use a WMAP cosmology in 
which (h, fl m ,n A ) =(0.7, 0.3, 0.7). 

2. OBSERVATIONS 

We used the Green Bank Telescope in 2007 to observe the 
z = 0.692 21 cm absorption spectrum against 3C 286. We 
observed all four Stokes parameters simultaneously using the 
digital FX Spectral Processor, which provides all the neces- 
sary self- and cross-products; here "FX" means that first it 
Fourier transforms the input signal and then multiplies the 
voltage spectra with appropriate phase shifts. This technique 
is described in detail by Heiles et al. (2001) and Heiles (2001). 

Our original paper reported a very statistically significant 
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detection of circular polarization (Stokes V) for the 21-cm 
line DLA system against 3C 286. As Zeeman splitting, it 
translates to a line-of-sight field strength of 84 ±9 ^G. Sub- 
sequently, after re-observing 3C 286 in January and March 
2009, we found that this detection is incorrect, basically be- 
cause what we reported as Stokes V — circular polarization — 
is really Stokes U — linear polarization. The discussion below 
explains the origin of this error, which is a matter of a mis- 
placed 90° of phase. 

2.1. Calibrating a Local Noise Source as a Secondary 
Calibrator Standard 

The sky electric field is sampled by probes in the telescope 
feed. In our observations with the GBT, these two native po- 
larizations are very close to being orthogonal linears. For 
convenience we denote the sampled voltage spectra (i.e., the 
Fourier transforms of the sampled voltages) as X and Y. Apart 
from small corrections for nonorthogonality and other cou- 
pling, the Stokes parameters are obtained as follows: 



I=XX+YY 
Q=XX-YY 
U = TRXXY) 
V = 21(XY) 



(1) 

(2) 
(3) 
(4) 



We tacitly assume that the voltage spectral products are time 
averages. Here 1Z and I mean the Real and Imaginary parts, 
respectively. To clarify these equations and their discussion, 
we neglect the small corrections for nonorthogonality; these 
corrections are embodied in the f eed's Mueller matrix and are 
determined as discussed in 32.31 

In a heterodyne radio receiver the signal from the sky is am- 
plified by mixing (multiplying) the sky signal by one or more 
locally-generated signals of known amplitude and phase. The 
resulting signal product, the gain G e , is a complex number in 
which the amplitude is the real part and the phase the imag- 
inary part. Between the receiver and the Spectral Processor 
there are a variety of electromagnetic components that alter 
the phase of the signal. Specifically, there are cables and opti- 
cal fiber, each of which adds phase proportional to 22£ where 
L is the length and A the wavelength of the signal in the partic- 
ular cable. Thus the gain G e is a complex number. Moreover, 
some of these electromagnetic components have their own po- 
larization properties that differ for X and Y, and as a result we 
have G e x and G e .y- For the complex portion of these gains, 
the only part of concern is the phase angle difference 



where 



tan< 



>e.Y : 



K(G e v) 



(5) 



(6) 



and a similar equation for X. 

These complex gains must be calibrated and their effects 
removed. The standard technique is to inject perfectly corre- 
lated noise at the feed from a noise source, commonly called 
the "cal". Like the sky signal, the correlated noise source in- 
jects voltages into the feed; unlike the sky signal, the cal's 
noise is perfectly correlated and corresponds to 100% polar- 
ization. We compare the cal's signal with that of a standard 



polarized calibration source of known flux to determine the 
equivalent antenna temperatures of the cal. We also compare 
the relative phase of the cal with that of the source, which 
relates the polarization of the cal to that of the sky. 

In essence, we determine the gain of the cal in X and Y 
relative to the sky signal; these gains are complex and change 
with frequency. We assume the cal's properties to remain con- 
stant with time. In subsequent observations, we use the cal's 
deflection to determine the complex electronic gains G e> x and 

2.2. Using the Cal to Calibrate the Complex Gains 

Consider the frequency dependence of equation [5] The in- 
dividual phase terms on the right-hand-side come from elec- 
tronic components and from cables. The phase delays of 
most electronic components change relatively slowly with fre- 
quency. However, for cables, the phase delay is and the 
phase difference in equation|5]is A</> = 2?r ^ L , where AL is the 
length difference between Y and X. In our experience, A(j) 
varies fairly rapidly with frequency. In particular, using the 
techniques described by Heiles (2001), for the current obser- 
vations we find 



dA(j> 
~dV 



;20 degMHz" 



(7) 



We have well-tested software that reliably performs this de- 
termination, and its correction (e.g. Heiles & Troland 2004). 
Unfortunately, however, this software had a bug, which we 
discovered on June 5, 2009 while we were analyzing both the 
2007 and 2009 data. If we fed the software an array of spectra 
to correct, it performed correctly. But if we fed the software 
just a single spectrum, it didn't apply any correction at all. Be- 
fore the current observations, we had never used the software 
on individual spectra. 

During our original 2007 observations this phase delay hap- 
pened to average about 91 degrees at band center with a slope 
of about 20 deg/MHz; the phase delay fluctuated by at most 
a few degrees for all of our spectra. For our bandwidth of 
625 kHz, the slope doesn't matter much. Thus, the phase 
difference of ~ 90 degrees effectively interchanges the real 
and imaginary parts of the cross-correlation product. In turn, 
this changes the derived Stokes V into Stokes U, and vice- 
versa. Fortunately, in our 2009 observations the phase delay 
was about 30 degrees, which alerted us to the error in the 2007 
data. Although the precise cause of the change in phase delay 
between 2007 and 2009 is unknown, we note that a change in 
cable length in either polarization over the entire signal path 
between the feed and correlator could result in such a change. 

2.3. The Mueller Matrix 

We corrected for polarization impurities in the system by 
using the Mueller matrix formulation given by Heiles et al. 
(2001). The Mueller matrix relates the measured auto- and 
cross-correlation products to the actual Stokes parameters. 
With this procedure, one observes a linearly polarized cali- 
brator over a range of parallactic angle PA to derive seven pa- 
rameters that describe the complex gains and coupling coef- 
ficients of the receiver system components. One then derives 
the Mueller matrix coefficients from algebraic combinations 
of these coefficients. 

In the present paper, we observed 3C 286. This source is, 
for most purposes at cm wavelengths, the "gold standard" 
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polarization calibrator. Specifically, the fractional polariza- 
tion and polarization position angle of 3C 286 have been sta- 
ble over four decades of observations across a wide range of 
frequencies above ~ 1400 MHz (Tabara & Inouye 1980; B. 
Gaensler 2009 priv. comm.). While these properties have 
not been established at lower frequencies, the dominance of 
the large-scale jet at these frequencies argue against any time 
variations of these quantities at 840 MHz. Hence we had the 
fortunate circumstance of observing a source that also hap- 
pens to be an excellent polarization standard. For about half 
of our observing days, we had enough PA coverage to derive 
the Mueller matrix from our observed data for that particular 
day; thus, we could use our observed astronomical data to also 
determine the Mueller matrix elements — a form of "self cali- 
bration". The derived corrections were closely the same from 
one day to another, although there were small differences. For 
any particular day's observations, we used the matrix that was 
closest in time to that day; for about half the days, this matrix 
was obtained on that very day. This gives us great confidence 
in our polarization calibration. 

Normally, one determines a single Mueller matrix for the 
observed band by averaging over the whole band. This leads 
to small calibration errors that change with frequency across 
the band, which in turn leads to small "baseline curvature" 
in the three polarization spectra. The latter are the fractional 
polarization p(v), polarization position angle x(v), and Stokes 
V parameter V(v); one normally fits a smooth baseline to these 
spectra and applies them as an ad-hoc correction. 

Here, we took a different approach. 3C 286 is such a 
strongly polarized source that, for any one day, we could 
easily determine the seven receiver system parameters on a 
channel-by-channel basis. We then fit their frequency de- 
pendencies by performing a minimum-absolute-residual-sum 
(MARS) 3rd-degree polynomial fifl We used these fit coef- 
ficients to derive the Mueller matrix independently for each 
channel. With this, the polarization calibration is very accu- 
rate for each channel independently, so that any features in the 
polarized spectra are real, characterizing the sky instead of the 
system. In particular, it is neither appropriate nor necessary to 
subtract off an ad-hoc baseline in any of the polarized spectra. 
It is these accurately-calibrated spectra that we use and plot in 
this paper. Even with this self-calibration of each channel in- 
dependently, there remain small residual inaccuracies, which 
result from applying Mueller matrices obtained on one day to 
data obtained on another day; this is probably responsible for 
the small spectral effects in the Stokes V spectrum in Figure 

m 

3. REANALYSIS 

First, we fixed the above-mentioned software bug. Then we 
reanalyzed the data. We used the same software as before and 
determined the Mueller matrices in the usual way (see Heiles 
et al. 2001), as we did before. We applied the Mueller matrix 
corrections and the phase difference correction of equation [5] 
as we did before. We did three things differently: 

1 . We derived the linear polarization by fitting the Paral- 
lactic angle dependence of Stokes U. Stokes U is a 

6 The MARS fit has the advantage of not responding to the 
spectral features caused by the absorption line. To our knowl- 
edge, this technique was invented by P.B. Stetson; see his website 

http://nedwww.ipac.caltech.edu/level5/Stetson/Stetson4.html HeileS haS addi- 
tional information on this technique; see Isf it_2008 . ps on his website 

http://astro.berkeley. edu /-he lies /handout s /handout s_num .html 



cross product and produces more reliable results, free 
of zero offsets, than Stokes Q, which is the difference 
between two large numbers. (As it happens, the re- 
sults from Stokes Q are indistinguishable from those 
of Stokes U). 

2. We reduced the off-line (i.e., the 3C 286 continuum) 
spectral channels to ensure that we obtain the correct 
results. That is, we did not "subtract off the continuum 
baseline". 3C 286 is the premier polarization calibrator 
in radio astronomy, and this ensures that our spectral 
line and continuum results correspond to what is known 
about 3C 286. 

The one wrinkle in this is that we have been unable to 
find a primary measurement of the polarization proper- 
ties of 3C 286 at frequencies as low as ours, 840 MHz. 
While 3C 286 exhibits essentially zero change with fre- 
quency of polarization position angle above 1 .4 GHz, 
we cannot be sure that the angle at 840 MHz isn't some- 
what different. We did not wish to address this question 
and therefore our reported polarization position angles 
have an arbitrary zero point. 

3. We edited the data much more carefully than before, 
which produced only minor improvements. 

4. THE PROPER DATA: STOKES PARAMETERS 

FigureQ]shows a four-panel plot of the data (black curves), 
now reduced properly. No "baseline corrections" have been 
subtracted or applied. The data consists of 12.6 hrs of on- 
source integration obtained in 2007. The orange curves are 
least squares fits to the data. The top panel shows Stokes 
I(v)/I (where / is the off-line continuum intensity, which 
is independent of velocity v), which has an absorption line 
with a fractional absorption of about 5%. The frequency 
centroid determined by the fit to the Stokes I spectrum is 
839.408348 ± 0.000046 MHz and its location is depicted 
by the vertical dot-dashed line. The corresponding redshift 
z=0.69215109±0.00000009 . The second pan el shows the per- 
cent polarization, p(v) = \J (U(v) 2 + Q( v) 2 ) / /, and the third 
panel shows the position angle x(v) (zero point is arbitrary). 

The linearly polarized p(v) line also appears in absorption 
(i.e. the percent polarization goes down in the line) and its 
fractional polarized absorption is close to that for Stokes /. Its 
line center differs from the Stokes / line center by 5.0 ± 0.25 
kHz. This 20<7 difference has high statistical significance. 
Taken together, these two results are consistent with the ab- 
sorbing gas covering both nonpolarized and polarized por- 
tions of the continuum source image, and the velocity of the 
part that covers the polarized portion differing from that cov- 
ering the unpolarized portion by 5.0 kHz (which is about 1.8 
km/sfl 

The third panel shows that the line changes the position an- 
gle by about 1.5 degrees at line center, which is consistent 

7 This shift is what looked like Zeeman splitting in our earlier result, where 
we thought that Stokes U was Stokes V. If the frequency shift for polarized 
intensity were reduced with the same software as we use to derive Zeeman 
splitting from Stokes V, then we would derive a field strength of 280 fj,G 
(instead of the original one in the paper of "only" 84 fj,G). This is because, 
when averaging linearly polarized Stokes parameters over hour angle, we 
properly rotate by the parallactic angle PA; in contrast, when we averaged 
our old Stokes V we didn't account for PA, so the amplitude was smeared 
down by improper weighting in PA). 
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Figure 1. Plots of Stokes parameters for 3C 286. From top to bottom, the panels are: the fractional absorption spectrum I(v)/I; fractional linear polarization 
p(v); position angle of linear polarization x( v ); ana fractional circular polarization V(v)/I. Fractional spectra are in percent and frequencies are bary centric in 
MHz offset from 839.408348 MHz. Black curves are data and orange curves are fits to the data. 
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with the absorbing gas covering only a fraction of the polar- 
ized source image and requires that the position angle of po- 
larization of the continuum changes across the source image. 

The fourth panel shows the percentage circular polarization 
V(v)/I. The absorption line exhibits a dip in fractional Stokes 
V of 0.03%. However, we believe this apparent circular polar- 
ization is meaningless. The line is 5% linearly polarized, and 
some of the total intensity and also the linear polarization will 
leak into Stokes V. The observed Stokes V intensity of the line 
is only about 1 % of its linearly -polarized intensity. Errors in 
the Stokes parameters at these levels are not unexpected, be- 
cause of the small day-to-day changes in the Mueller matrix 
elements. 

5. CLOUD STRUCTURE PRODUCING THE 
ABSORPTION FEATURE 

We now show that the GBT polarization spectra provide 
new and unique information about the kinematic and spatial 
structure of the gas that gives rise to 21 cm absorption to- 
ward 3C 286. In particular we note the following model- 
independent facts: (1) The difference between the velocity 
centroids of the Stokes / and fractional polarization spectra is 
inconsistent with a single-cloud origin for the absorption. (2) 
The Gaussian shape of the absorption profile is plausibly ex- 
plained through the cental limit theorem by the superposition 
of many small clouds. (3) The position-angle spectrum rein- 
forces the conclusions drawn from the fractional polarization, 
but requires detailed modeling as described below. Besides 
studying the spatial and kinematic structure of the absorbing 
gas, we also investigate its two-phase structure; i.e., whether 
the gas is a cold neutral medium (CNM) with temperature T 
~ 100 K, a warm neutral medium (WNM) with T ~ 8000 
K (Wolfire et al. 1995), a thermally unstable phase with tem- 
perature between these two extremes, or some combination of 
any of these. Furthermore, we consider clues to the identity 
of the galaxy hosting the absorbing gas. 

5.1. Core-jet structure of the radio source 

The 21 cm absorption likely forms against a source with 
the core-jet configuration shown in Fig. [2] (Wilkinson et al. 
1979). This VLBI map, obtained at 609 MHz, indicates that 
3C 286 is an asymmetric source consisting of a compact core 
and an extended jet. The core, which covers a solid angle, 
^core~10x5 milli-arcsec 2 (mas 2 ) with major-axis position an- 
gle equals 42°, contributes a flux density 5'„(609)=3.5 Jy, 
while in the case of the jet Qj et =40x 15 mas 2 , PA.=42°, and 
5'„(609)=14 Jy. Since a similar core-jet structure is detected 
at 329 MHz, 1667 MHz (Simon et al. 1980), 5 GHz (Jiang et 
al. 1996; Cotton et al. 1997), and 15 GHz (Kellermann et al. 
1997, priv. comm.), it is safe to assume that such a structure 
is present at the absorption frequency of 839.4 MHz. Adopt- 
ing the spectral indices measured by Simon et al. (1980) of 
cw = -0.29±0. 15 and a jet = -0.55, we find that S™ re «3.2 Jy 

and sjfwl 1.7 Jy at v k, 840 MHz. 

While the polarization structure of 3C 286 has not been 
detected at low frequencies, VLBI measurements of all four 
Stokes parameters have been obtained at 5 GHz (Jiang et al. 
1996; Cotton et al. 1997; Cotton 2010 [priv. comm.]). These 
data show that unlike most quasar VLBI sources, the frac- 
tional polarizations, i.e., fraction of total surface brightness in 
y/U 2 + Q 2 , where U and Q are the Stokes parameters, of the 
core and jet are comparable. More specifically, the ratio of jet 
to core polarized surface brightnesses is about 0.4. 




N 



Figure 2. VLBI map of 3C 286 at 609 MHz (Wilkinson et al. 1979). Con- 
tinuum structure exhibits an unresolved core source to the north east and an 
extended jet source to the south west. Orthogonal straight lines are 0.02 arc- 
sec in length corresponding to 142 pc at the absorber. 

Nevertheless the ratio of polarized flux densities P' et /p core ss 
5 owing to the much larger solid angle subtended by the jet, 
where the polarized flux density P = \/U 2 + Q 2 . Because of 
the steeper spectral index of the jet, it is reasonable to as- 
sume that this ratio increases with decreasing frequency. As 
a result, we shall ignore the core source when computing ab- 
sorption in Stokes U and Q parameters, or equivalently in the 
spectrum of polarized flux density P (note, P = pi) and po- 
larization position angle x = 0.5 x arctan(C/ /Q). We check the 
self-consistency of this assumption in § 7. 

5.2. Two Cloud Model 

We now describe the procedures used to model the Stokes 
parameter spectra. 

5.2.1. Model Description 

To account for the absorption spectrum we first adopt a sim- 
plified model, illustrated in Fig. [5] comprised of one uniform 
H I 'cloud' toward the jet (cloud 1) and another uniform H 




N 



Figure 3. Foreground clouds (grey) superposed on VLBI map from Fig. [5] 
Cloud 2 at NE completely covers core source, while cloud 1 toward SW par- 
tially covers resolved jet source. Cloud 1 is shaded darker to illustrate its 
larger 21 cm optical depth. Light blue shading indicates that region emitting 
polarized radiation is the jet. Thick red bars depict direction of E vector po- 
larization. Orthogonal straight lines are 0.02 arcsec in length corresponding 
to 142 pc at the absorber. 
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I 'cloud' toward the core (cloud 2). A single cloud model is 
ruled out by past VLBI observations of the 839.4 MHz ab- 
sorption feature that indicated an inhomogeneous structure of 
the absorbing gas. The velocity difference between the max- 
imum depths of the phase-shift and fringe amplitude spec- 
tra detected by a two-element VLBI experiment (Wolfe et al. 
1976) indicates a velocity difference between the H I towards 
the core and jet sources. As a result the relative change in 
Stokes / parameter at velocity v is given by 



A/(v)// = £)q/ i 



exp - r,(v) - 1 



(8) 



where the number of clouds n c =2, C\ and C2 are the area cov- 
ering factors of the gas toward the jet and core, f\ and fz are 
the fractions of Stokes / flux density in the jet and the core 
sources, and t\(v) and T2(v) are the 21 cm optical depths av- 
eraged across the areas subtended by jet and core at the ab- 
sorber. Note, / in the denominator of eq. (4) is the continuum 
Stokes / parameter extrapolated across the 21 cm absorption 
line. 

In order to fit the model to the data we made the following 
assumptions (see Fig.0: First, because most of the polarized 
flux density arises from the jet, we assume the absorption 
spectrum of polarized flux density, P(v), arises only in 'cloud' 
1 . Since the line center of the P(v) spectrum is shifted by + 
5.0±0.25 kHz (-1.8±0.1 km s" 1 ) relative to the line center 
of the Stokes /(v) spectrum, we assume that the velocity cen- 
troid of cloud 1 is lower than that of cloud 2. Second, the shift 
in polarization position angle with frequency , Ax(v)=x(v) - 
Xcont (where x cont is the polarization position angle averaged 
over the continuum source), indicates that 'cloud' 1 partially 
covers source 1, i.e., C\ < 1, and the intrinsic value of x varies 
across the source, as is observed in VLBI experiments at high 
frequencies (Jiang et al. 1996; Cotton et al. 1997). With a 
covering factor C\ < 1, the position angle of polarized flux 
behind the cloud is weighted less in the line than in the con- 
tinuum, thereby causing a shift of position angle in the line. 
Third, because of the smaller dimensions of the core source, 
we assume Cz=l. Fourth, we assume Gaussian velocity distri- 
butions for each cloud, in which case the optical depth of the 
; th cloud at velocity v is given by 



Ti(v) = T 0i ,xexp 



((v- Vi )/V2a v ,ij 



(9) 



where v,- and and er,,,- are the velocity centroid and velocity 
dispersion of the i' h cloud. In the case of 21 cm absorption 
the central optical depth 



(10) 



T ,i = N m ,/ U^xlO'V,, t T„ 



where T s j is the spin temperature and ex,.., is in units of km 
s -1 . Fifth, because the polarized radiation is assumed to be 
emitted only by the jet source, the relative change in fractional 
polarization is given by 



AP(v)/P = d 



exp -ti(v) -1 



(11) 



Because the source is weakly polarized, the flux from the jet is dominated 
by Stokes / with smaller contributions from Stokes U and V comprising the 
polarized contribution. 



where we again extrapolate across the line to obtain P in the 
denominator of the last equation. The geometry of the model 
is illustrated in Fig. [3] 

5.2.2. Model Fits To The Data 

We adopted an incremental approach by first modeling the 
Al /I spectrum alone. When a successful fit was found we 
then fitted Al/I together with either the AP/P or Ax spec- 
trum simultaneously. If a successful fit were obtained, we 
followed up by fitting Al/I, AP/P, and A\ spectra simulta- 
neously. The 1-cr errors, which were derived by computing 
the standard deviations from portions of the spectra displaced 
from the absorption features, are as follows: cta///=0. 00050, 
<j ap/ p=0.0034, and cr Ax =0.042°. 

We obtained the relevant parameters by fitting the model 
to the data with the Levenberg-Marquardt method that makes 
use of the non-linear chi square minimization routine, mrqmin 
(Press et al. 1996). For the Al /I spectrum we fixed the 
area covering fractions Ci=0.5 and C^= 1 .0 as discussed above, 
but allowed the other parameters to float after making initial 
guesses that were guided by the discussion in § 5.2.1. Specifi- 
cally we solved for Nms /T s j 0, v,-, er v .i, and /;. A successful fit 
was found and the output parameters were then used as inputs 
for fitting the Al /I and AP/P spectra simultaneously. The 
results are shown in Fig. |4] 

While the fit to the Al /I spectrum appears to be quite accu- 
rate, the fit including the AP/P spectrum has a chi square per 
degree of freedom, x 1 /v=YL.(s, which is unacceptable. The 
reasons for the failure of the model are straightforward. Both 
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Figure 4. Stokes spectra of 3C 286. Solid curves are two-cloud fits and data 
given by histograms. Top panel shows A/// spectrum. Spectrum with dashed 
curve showing solution for cloud 1 and dot-dashed curve for cloud 2. Bottom 
panel shows AP/P. 

9 Since we make no a priori assumptions about the value of Nm i m either 
cloud, we cannot solve for Njjh nor T s ; separately. 
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the depth of the model AP/P absorption feature and the loca- 
tion of its frequency centroid are in significant disagreement 
with the data. Altering the input parameters to get a better 
fit to the AP/P spectrum does not solve the problem, be- 
cause such changes also alter the Al/I spectrum, which result 
in significant disagreement between the model and observed 
Al/I spectra. The dilemma is that while cloud 1 produces 
both the AP/P spectrum and most of the Al/I spectrum, the 
difference between the velocity centroids of these spectra are 
too large to be modeled with a single velocity component. Be- 
cause of a similar difference between the frequency centroids 
of the Ax and Al/I spectra, it is clear that the polarized radi- 
ation is incident on a third velocity component. As a result we 
abandon the two-cloud model in favor of a three-component 
model to which we now turn. 

5.3. Three Component Model 

For this model we again ignore polarized radiation from the 
core and assume that it is emitted only by the jet. Again 
we assume that most of the flux from the jet is in Stokes /. 
But we now assume that part of the jet has low or vanishing 
polarization and that 'cloud' 1 covers both this unpolarized 
portion of the jet and that part of the jet which emits polar- 
ized radiation (the faint blue section of the source in Fig. [5]): 
the covering factor for the latter is C3. We further assume the 
presence of a velocity gradient in 'cloud' 1 such that the ve- 
locity centroid for the Stokes / spectrum is again Vi, but for 
the AP/P spectrum is now V3 (see Fig.[5j. As a result we now 
have 



AP(v)/P = C 3 



exp -r 3 (v) -1 



(12) 



where tj,(v) is the 21 cm optical depth of the gas in cloud 1 
toward the polarized portion of the jet source: the velocity 




v X 



Figure 5. Same as Fig.f5]except that cloud 1 covers both unpolarized (white) 
and polarized (blue) portions of the jet. Gas in cloud 1 toward polarized 
portion is called component 3, which has a velocity centroid lower than that 
of the entire cloud 1 . Note that intrinsic polarization position angle \ (where 
X is defined to increase positively in counter-clockwise direction from due 
North) decreases along the source in the SW direction. x = Xb for the area of 
the polarized jet behind cloud 1 (i.e., behind component 3) and \ = \ a for the 
area of the jet not behind cloud 1. Orthogonal straight lines are 0.02 arcsec 
in length corresponding to 142 pc at the absorber. 

10 Although the core is observed to be polarized at higher frequencies [e.g. 
Jiang et al. 1996], in §5.1 we argued that the flux density of polarized radia- 
tion mainly comes from the jet. 



centroid of this gas is V3. Therefore, this model has velocity 
components centered on vi , V2, and 1/3. 

Since we shall also consider fits to the Ax spectrum in this 
case, we develop a model for the intrinsic change in polar- 
ization position angle across the jet source. We assume that 
X = Xa for the polarized flux density emitted by the fraction of 
the jet source that is not incident on the cloud and x = Xb for 
the fraction of emitted polarized flux density that is incident 
on the cloud. Noting that the net flux densities in Stokes Q 
and U parameters are given by Q = Q a + Qb and U = U a + Uf,, 
and x=0.5 xarctan(/7 / Q), we find that 



X(v)= ^arctan 



rC 3 tan(2xz,)exp(-T3(v)) + ( 1 - C 3 )tan(2x fl ) 



rC 3 exp(-r 3 (v)) + 1 - C 3 

(13) 

Here r=< Xq >/,/< Iq > a , where < Iq >„ and < Iq >& 
are the surface brightnesses in Stokes Q averaged over re- 
gions a and b respectively, and we make use of the relation 
Qb/Qa = rCi/f/L—C-}). The observed value of x averaged over 
the unattenuated continuum source, Xcont, is obtained from the 
last equation by setting T3=0. The difference between x in the 
line and Xcont is denoted by Ax(v)=x-Xcont- 

We fitted the three-component model parameters to the 
spectra as follows. First, observing that the parameters char- 
acterizing the Stokes Al/I spectrum are independent of those 
characterizing the AP/P and Ax spectra, we did not fit the 
three Stokes spectra simultaneously. Rather, we fitted the 
Al/I spectrum alone, and then the AP/P and Ax spectra si- 
multaneously. Given the reasonable fit to the Al /I spectrum 
discussed in § 5.2, we adopt the parameters found for clouds 
1 and 2 in that fit here. Second, in the case of the polarized 
spectra we need to find the gas parameters, A^Hij/Ti^, V3, av.3, 
and C3 as well as the source parameters r, Xa> and Xb- For 
the source parameters, we are guided by single-dish polariza- 
tion measurements above 1.5 GHz that show the continuum 
polarization position angle for the entire source Xcont=33±5° 
(Tabara and Inoue 1980). The VLBI polarization data at 5 
GHZ (Jiang et al. 1996; Cotton et al. 1997; Cotton 2010 
priv. comm.) further show that along the jet x decreases with 
increasing distance from the core source (i.e., in the SW di- 
rection) by as much as w 40° which is consistent with the 
orientation of the polarization bars in Fig [5] We found that 
the Ax(v) line spectral feature is more naturally explained by 
a large central optical depth to 3 combined with a low cover- 
ing factor C3 rather than vice versa. By contrast, C3 and T0.3 
are degenerate in the case of the AP/P spectrum (see eq. (8)). 

We used mrqmin to determine the cloud 1 and cloud 2 pa- 
rameters by fits to the Al /I spectrum. As discussed above, 
the parameters of the third velocity component were found by 
separately fitting the AP/P and Ax spectra together. This 
is possible since the polarized spectra are independent of the 
cloud 1 and 2 parameters, while the Stokes Al /I spectrum is 
independent of the component 3 parameters. After trials with 
several values of C3 we fixed C3=0.2 at the outset. Larger val- 
ues of C3 must be accompanied by smaller values in 70,3 in 
order to maintain the observed value of AP/P (see eq. 8), 
but an increase in C3 generates unacceptably large values of 
Xb~Xa in order to retain the maximum dip of Ax~ -1.5°(see 
eq. 9). On the other hand smaller values of C3 imply larger 
values of the ratio < Iq >b / <Iq > a , i.e. larger variations of 
polarized surface brightness across the source, than observed 
at 5 GHz (Jiang et al. 1996; Cotton et al. 1997). We consid- 
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Cloud 



Tqi 



Nm.,IT SJ 
crrf J K H 



Gas Properties 



Source Properties 



km s 



km 



-l 



fi 



Xa 
deg. 



Xb 
deg. 



0.128±0.005 
0.057±0.005 
0.280±0.004 



(2.03±o.oi)xio jS 

(7.94±0.10)xl0 17 
(4.40±0.02)xl0 ls 



-0.82±0.02 3.48±0.02 0.5 
+2.37±0.07 3.05±0.07 1.0 
-1.44±0.04 3.42±0.04 0.2 6 



0.72 
0.28 

>0.40 c 2.70±0.9 



43.6±0.9 



Table 1 

Properties of Clouds and Source Causing the Absorption Feature 



(a) Velocity v=0 corresponds to 2=0.69215109, shown as dotted vertical lines in Figs. 1,4, and 6. 

(b) Fraction of polarized flux density incident on "cloud 3". 

(c) Fractional area of jet source emitting polarized radiation, assuming r=l, 



ered letting C3 be one of the floating variables, but decided to 
let it remain fixed, because the spatial variation of \ across the 
jet is at best a qualitative rather than a robust constraint. For 
these reasons we fixed r=l at the outset (where r is the ratio 
of average Stokes Q surface brightnesses in region b relative 
to a), but let the remaining parameters float. 

The results are shown in Fig. [6] which plots model and ob- 
served spectra. The output parameters are listed in Table 1 
along with 1-er errors. These were obtained from Ax 2 confi- 
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Figure 6. Stokes spectra of 3C 286. Solid blue curves are three- 
component fits and data are magenta histograms. Top panel shows 
AI/I spectrum with cloud 1 and 2 solutions as dashed and dot- 
dashed black curves. Middle and bottom panels show AP/P and 
Ax spectra. Model parameters are in Table 1. Vertical dotted 
line at offsets =0 corresponds to centroid of Stokes / feature at 
z=0.69215109±0.0000000109. Dashed vertical line in top panel at 
v=4.21±0. 1 1 km s~ shows velocity predicted by the optical redshift 
atz=0.69217485±0.00000058 (Wolfe et al. 2008). 



dence ellipses containing 68% of the distributed data by pro- 
jecting the ellipses onto each parameter axis (see Press et al. 
1996). We did this for all possible pairs of the parameters in 
Table 1 and chose the largest value of the projection as the 
l-cr error. While the fit to the Al /I spectrum appears to be 
quite good, the chisquare per degree of freedom x 2 /V=1.73, 
which corresponds to a probability of salO -8 that \ 2 exceeds 
294 for 169 degrees of freedom. One reason for the elevated 
X 2 / v is that our model is undoubtedly inadequate given the s» 
100:1 signal-to-noise ratio of the Al /I data. In principle we 
could improve the fit by adding more velocity components. 
But this would be an unconstrained ad hoc procedure, given 
our limited knowledge of the source-gas configuration, and 
thus we decided against it. The main reason for the high value 
of x 2 l v is that the largest contribution to x 2 comes from the 
relatively large scatter of the data around the model between 
839.385 MHz and 839.397 MHz (i.e., between Av=- 0.023 
and - 0.01 1 MHz in Fig. 6). As the cause of this effect is cur- 
rently uncertain, we attribute it to some unknown systematic 
error. 

The fits to the polarized spectra result in x 2 /V=1.54. This 
corresponds to a probability of sa 3xl0~ 10 that x 2 exceeds 
536 for 348 degrees of freedom. The reasons for the high 
value of x 2 l v m this case stems from the small, m 3 kHz, 
difference between the velocity centroids of the AP/P and 
Ax spectra. In contrast to the velocity difference between the 
centroids of the polarized and Stokes / spectra, this effect is 
difficult to understand. In any plausible model the velocity 
centroids of the AP/P and Ax spectra coincide with the cen- 
troid of the optical-depth velocity distribution (see eqs. 4, 8 
and 9). We tested the possibility that the error might be due to 
a constant frequency offset between the AP/P and A% pro- 
files by shifting the frequency scale of the AP/P spectrum. 
We found that x 1 l v decreased to a minimum of 1.48 when 
the shift was one frequency pixel (i.e., 1.2 kHz). However, 
the improvement of the fit was insufficient for the error to be a 
frequency offset. More likely this is also an unknown system- 
atic error. Adding another velocity component doesn't help 
since though it might improve the fit to the AP/P spectrum 
at 839.405 MHz, it also worsens the fit to the AI/I spectrum. 
Despite these problems, the model provides an adequate fit to 
the data, given that our goal is to determine a general picture 
of the cloud structure rather than obtain a precise evaluation 
of cloud sizes, locations, etc. As a result we shall adopt the 
three-component model in further discussions. 

6. DISCUSSION 

In this section we discuss implications of the models used 
to describe the 21 cm absorption in 3C 286. 

6.1. Cloud and Galaxy Properties 
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In § 5.3 we discussed the three-component model used to 
model the AT/7, AP/P, and Ax absorption spectra. As illus- 
trated in Fig. [5] the three-component model is actually based 
on two clouds. One cloud (cloud 2) completely covers the 
core source. The other cloud (cloud 1) partially covers the 
jet source, and contains a velocity gradient such that the part 
of the cloud toward the polarized portion of the jet (compo- 
nent 3) causes the velocity centroid of the polarized spectra 
to be offset in velocity from the unpolarized A7/7 spectrum. 
We simplified the model by letting component 3 be a uniform 
sub-region of cloud 1 with a velocity centroid, velocity dis- 
persion, and 21 cm optical depth that differ from their values 
in the parent cloud. Table 1 indicates that the polarization 
position angle for the part of the source absorbed by the fore- 
ground cloud is Xi=43.6°. Because the ratio of the Stokes 
parameters, U / (2=tan(2x), the predicted ratio U / (2=20.4. By 
contrast, U /Q=0.046 toward the part of the jet that bypasses 
the absorbing gas where the model predicts x«=2.67°. This 
is consistent with the independent Stokes U and Q absorp- 
tion spectra, which show a clear detection of the absorption 
feature in Stokes U but not in Stokes Q. The absorption fea- 
ture, which is detected at a signal-to-noise ratio of about 20: 1 
in Stokes U, is predicted to be below the 1-er noise level in 
Stokes Q. The model also predicts Xco«r=38° for the position 
angle integrated over the entire source, which is consistent 
with single-dish measurements of 33±5° at 1.5 GHz (Tabara 
& Inoue 1980). Our model predicts that the intrinsic polar- 
ization position angle of the jet decreases by \a - Xb= -40.9" 
in the SW direction (see Fig . [5j> ■ Polarization maps at 5 GHZ 
(Jiang et al. 1996; Cotton et al. 1997) do show shifts by 
about this amount in the sense of decreasing x along the jet 
axis in the SW direction. However, the value of x at the ex- 
treme edge of the jet appears to be greater than our prediction 
of x a =2.67 . These predictions could be checked with VLBI 
polarization measurements at lower standard frequencies such 
as 609 MHz. 

Next we turn to the properties of the absorbing gas. The 
large projected areas subtended by the sources at the ab- 
sorber (A core ~70x35 pc 2 for the core and Aj et w280x60 pc 2 
for the jet) likely indicate that the gas causing 21 cm absorp- 
tion in 3C 286 is comprised of many interstellar clouds. In 
the Galaxy a beam with area Aj et directed perpendicular to 
the plane of the disk would subtend A/"cnm~1000 CNM (i.e., 
cold neutral-medium [Wolfire et al. 1995]) clouds. This is 
because A/cNM=«CA , MxAj e tx2// where the space density of 

CNM clouds n C NM ~ 3x 10" 4 pc" 3 (McKee & Ostriker 1977) 
and the disk half-thickness H ss 100 pc. This value of A/cnm 
is a lower limit, since A/'cnm=«ca'm x Ax277/cos(0 for a disk 
with inclination angle i. For the core source A/"cnm>60. As 
a result the Gaussian shape of the absorption profile is natu- 
rally explained as the optical-depth weighted sum of multiple 
Gaussians (i.e., the central limit theorem). 

Note that the central velocities of the Gaussians are likely 
superposed on large-scale velocity gradients present both in 
'cloud' 1 toward the jet and 'cloud' 2 toward the core. While 
such a gradient in 'cloud' 1 is required to explain the velocity 
shifts between the stokes 7 and polarized spectra, as explained 
above, the necessity for a gradient in 'cloud' 2 stems from the 
approximate agreement between the optical redshift and red- 
shift of 'cloud' 2 illustrated in Fig. 6. This implies that the 
optical continuum source is, as expected, physically associ- 
ated with the compact radio core. The small, but significant, 
difference between these redshifts further suggests that since 



the optical beam size is small compared to the dimensions of 
'cloud' 2, the optical beam samples only a limited portion of 
the velocities spanned by the gradient in 'cloud' 2. Therefore, 
although the UV resonance lines form in gas within 'cloud' 
2, the velocity centroids formed by averaging over the optical 
and wider radio beams need not be equal. 

To determine whether the absorbing gas is CNM, WNM, or 
something else we need to determine its kinetic temperature, 
Tk- Assuming the foreground gas covers the entire source with 
a uniform cloud characterized by <r v = 3.75 km s" 1 (Wolfe et 
al. 2008), one finds an upper limit of r™ x =1690 K. If we 
further assume that Nm equals the UV-determined value of 
1.77xl0 21 cm" 2 (e.g. Wolfe & Davis 1978; Kanekar 2003a), 
then eq. (6) implies 7^=1035 K, since the central 21 cm optical 
depth T =0.10 (Wolfe et al. 2008). Because T s <T k < T™* 
for diffuse warm gas (Liszt 2001), these results would indicate 
the gas is neither standard CNM where Tjss 80 K nor standard 
WNM where T k « 8000 K (Wolfire et al. 1995), but rather 
resembles the thermally unstable phase found by Heiles & 
Troland (2003) in the Galaxy ISM. 

On the other hand the frequency-dependent change in po- 
larization position angle across the absorption feature (panel 
3 in Fig. |5]l provides strong evidence against the assumption 
of a single uniform cloud (see § 5.4). From the values of <r v ; 
in Table 1 we find that r™ x =1450 K, 1115 K, and 1400 K 
for gas in velocity components 1, 2, and 3. If we again as- 
sume that A^hi,;=1.77x 10 21 cm" 2 in each case, the values of 
Nm,ilT s ,i in Table 1 indicate that 7; u =892 K 7^=2230 K, and 
r s .3=394 K. However, since the values of Nm,i/T S i in Table 1 
are spatial averages over dimension much larger than the ss 1 
It. yr. size of the UV continuum, it is unlikely that the UV 
determined value of Nm applies to the gas in each compo- 
nent, nor that each component has the same value of Nhi- It 
is equally plausible to assume that A^hi averaged across cloud 
1 is a factor of 2 lower than the UV determined value, and 
C\ equals 0.25 rather than 0.5. In that case r, i=223 K. Sim- 
ilarly if Nm3=Nm.\, then 71.3=197 K. Since these are physi- 
cally plausible temperatures for CNM gas with the low metal 
abundances inferred for this absorber (Wolfe et al. 2008), it is 
reasonable to assume they represent kinetic temperatures. In- 
terestingly, the assumption A r Hi.2=l-77x 10 21 cm" 2 results in 
the unreasonable prediction that T s .2 is higher than 7^™": we 

take this as direct evidence that Nm.2 is lower than 1 .77 x 10 21 
cm" 2 and/or C% < 1, thereby implying that cloud 2 could also 
be CNM. Of course, we cannot entirely rule out the possibil- 
ity that the gas is in a thermally unstable phase. But the above 
arguments in addition to the results from the VLBI line ex- 
periment, which indicates individual components toward the 
core and jet with kinetic temperatures T k < 500 K (Wolfe et 
al. 1976), make a compelling case that the gas causing 21 cm 
absorption in 3C 286 is mainly CNM. Although this conclu- 
sion is at odds with the anti-correlation between T s and [M/H] 
found by Kanekar et al. (2009b), the evidence provided by 
our polarization measurements suggests either that physical 
conditions in this absorber deviate from typical conditions in 
DLAs or that most values of T s deduced from UV-determined 
values of A^hi are poor indicators of T k . 

Identification of the galaxy hosting the absorbing gas is dif- 
ficult to ascertain. Le Brun et al. (1997) used HST images 
to tentatively identify an object (2c in their nomenclature) 
as the leading candidate for the host. Although its absolute 
magnitude Mb = -20.1 (in our cosmology) suggests a mas- 
sive galaxy, the metallicity [M/H]=-1.6 of the absorbing gas 
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(Wolfe et al. 2008) is at the low end for this redshift, which 
suggests a low-mass galaxy (Tremonti et al. 2004). Although 
the area subtended by 3C 286 at the absorber would encom- 
pass over 1000 CNM clouds, the 3.75 km s" 1 velocity dis- 
persion of the absorption feature more closely resembles that 
of a single cloud rather than an ensemble of clouds. But this 
might be explained by a reduced supernovae input of mechan- 
ical energy into gas comprised of many clouds. This is consis- 
tent with the low metallicity and absence of C II* absorption 
both of which indicate a lower than normal SFR (Wolfe et al. 
2008). As a result, the evidence accumulated so far suggests 
that the gas causing 21 cm absorption in 3C 286 is embedded 
in a massive galaxy with a low SFR. 

6.2. Limits on Variations of Physical Constants 

The difference between the optical redshift, z opt , and radio 
redshift, z, a dio, places limits on variations of physical con- 
stants between the absorption epoch and the present, which 
is given by 

A\n(a 2 g p m/M) = |(z pt-Z ra dio)/(l +&adio)|, (14) 

Here a is the fine structure constant, g p is the gyromag- 
netic ration of the proton, and m/M is the electron-to- 
proton mass ratio (Wolfe et al. 1976). Because the ra- 
dio redshift z radio =0.69215 109±0.0000000109 and the optical 
redshift z opt =0.69217485±0.00000058, A\n(a 2 g p m/M)=A.2 
km s" 1 ±0.10 km s" 1 (in velocity units) where the error is 
dominated by the uncertainty in z opt . However, this estimate 
ignores the much larger systematic error determined from 
night-to-night changes in wavelength calibration on HIRES, 
which can be as large as 2 km s" 1 (Kanekar et al. 2010). 
Therefore, adopting the 4.2 km s" 1 difference as a conser- 
vative 2cr upper limit, we find Aln(a 2 g p m / M) < 1.4xl0~ 5 
in the redshift interval z=[0, 0.692]. On the other hand our 
conclusion that the optical source is located in the compact 
core source implies that the relevant radio redshift is given by 
cloud 2, which is shifted by 1.8 km s" 1 from z opt (see Table 
1). In that case Aln(a 2 g p m/M) < 0.6xl0" 5 . Both of these 
limits are comparable to results from recent measurements by 
Kanekar et al. (2010; see also Tzanavaris et al. 2007) who 
used redshifts deduced from C I absorption to determine the 
optical redshift, since C I is a more accurate tracer of the CNM 
gas that gives rise to 21 cm absorption than the Si II, Fe II, Zn 
II, and Cr II resonance transitions used to determine the opti- 
cal redshift of the 21 cm absorber toward 3C 286 (Wolfe et al. 
2008). In principle the present limit could be improved with a 
future space based detection of C I absorption lines in this ab- 
sorber, and a reduction of the systematic errors in wavelength 
calibration of the HIRES spectrograph. 

7. SUMMARY AND CONCLUSIONS 

Reanalysis of the polarization spectra of the 21 cm absorp- 
tion line detected at z=0.692 toward 3C 286 leads to the fol- 
lowing conclusions: 

(1) The detection of an 84 /iG magnetic field inferred from 
Zeeman splitting by Wolfe et al. (2008) is erroneous. An in- 
sidious software error failed to remove a spurious phase delay 
of pa 90° between orthogonal linearly polarized signals from 
the GBT antenna. This phase difference changed the derived 
Stokes V parameter into Stokes U and vice versa. As a result, 
what we interpreted as the spectral signature of Zeeman split- 
ting in Stokes V was a shift between the velocity centroids of 



the Stokes / and Stokes U absorption spectra. 

(2) A proper reanalysis of the GBT spectra is presented in 
Fig. 1. The Stokes / spectrum again exhibits the 21 cm ab- 
sorption feature at v = 839.408348 MHz (z=0.6 921 1510 9). 
The spectrum of fractional polarization, p=-\/ Q 2 + U 2 /I, 
where Q and U are the Stokes parameters, shows that the 
frequency centroid of the absorption feature is offset by 
5.0±0.25 kHz from the frequency centroid of the Stokes / 
spectrum. The polarization position-angle spectrum exhibits 
a frequency centroid that is also offset from the frequency cen- 
troid of the Stokes / spectrum, but by an amount that differs 
slightly from the centroid of the fractional polarization spec- 
trum. The Stokes V spectrum exhibits an absorption feature 
at 839.41 MHz, but this is most likely an artifact due to leak- 
age from the other Stokes parameters. The latter imply a 3-rr 
upper limit of 17 /zG on the line-of-sight component of the 
magnetic field in the absorbing gas. 

(3) We modeled the absorption feature with a core-jet ra- 
dio structure of 3C 286 behind three velocity components: in 
§ 5.2.2 we describe why two velocity components will not 
work. In the three-component model, cloud 2 completely cov- 
ers the core source and cloud 1 partially covers the jet with a 
sub-region 3 that partially covers the polarized portion of the 
jet. The radio beam through cloud 1 has a covering factor 
C\ for the Stokes / radiation, but a smaller covering factor C3 
for the polarized radiation (see Fig. |5). Since we assume the 
presence of a velocity gradient in this cloud the velocity cen- 
troids in Stokes / and the polarized radiation in the jet will dif- 
fer. The predictions of this model are in reasonable agreement 
with the data. Comparison between the model and observed 
spectra is shown in Fig. [6] and the physical parameters of the 
model are given in Table 1 . 

(4) In §7 we discuss implications of these models. We show 
that the change in polarization position angle predicted by the 
model is consistent with polarization measurements made at 
higher frequencies, but with some important differences that 
may be tested with future VLBI observations. We show that 
the area of the jet projected onto the absorbing gas would en- 
compass more than 1000 standard CNM type clouds found in 
the Galaxy. This provides a natural explanation for the Gaus- 
sian shape of the absorption feature through the central limit 
theorem. We argue that the high spin temperatures deduced 
from the single uniform cloud model are unlikely to be cor- 
rect due to evidence for an inhomogeneous distribution of the 
foreground gas. While we cannot rule out the possibility that 
we have detected warm, thermally unstable gas, we have pre- 
sented independent arguments favoring the CNM hypothesis. 
The evidence accumulated so far also indicates the absorbing 
gas is in a massive galaxy with a lower than average SFR. 

(5) A comparison between the radio and optical redshifts 
of the 21 cm absorber sets a conservative 2a upper limit on 
variations on the product of three physical constants given 
by Aln(a 2 g p m/M) < 1.4xl0" 5 within the redshift interval 
z=[0,0.692]. Based on our conclusion that the optical absorp- 
tion occurs in cloud 2, this limit reduces to Aln( a 2 g p m/M) < 
0.6 xlO" 5 . 

To conclude, our measurement of 21 cm absorption in all 
four Stokes parameters has provided new information about 
the foreground gas at z=0.692 toward 3C 286. The difference 
between the velocity centroids of the absorption feature in un- 
polarized radiation on the one hand and in fractional polariza- 
tion and polarization position angle on the other hand demon- 
strates evidence for spatial variations of Nm/T s on scales less 
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than ~ 100 pc. Our observations demonstrate that simple 
models of uniform clouds covering the entire radio source 
structure are incorrect. This is clearly illustrated by our anal- 
ysis of cloud 2 in front of the compact core source: applica- 
tion of the UV-determined value of Nui leads to a value of 
the spin temperature significantly higher than the upper limit 
to the kinetic temperature set by the velocity dispersion of 
the gas, which is physically implausible. Instead the data are 
more consistent with a model in which each of the three veloc- 
ity components is comprised of large numbers of small CNM 
clouds that produce the smooth Gaussian velocity profile of 



the absorption feature. 
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